Theory of the [111] magnetization plateau in spin ice 
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The application of a magnetic field along the [111] direction in the spin ice compounds leads to two 
magnetization plateaux, in the first of which the ground state entropy is reduced but still remains 
extensive. We observe that under reasonable assumptions, the remaining degrees of freedom in the 
low field plateau live on decoupled kagome planes, and can be mapped to hard core dimers on a 
honeycomb lattice. The resulting two dimensional state is critical, and we have obtained its residual 
entropy - in good agreement with a recent experiments - the equal time spin correlations as well as 
a theory for the dynamical spin correlations. Small tilts of the field are predicted to lead a vanishing 
of the entropy and the termination of the critical phase by a Kasteleyn transition characterized 
by highly anisotropic scaling. We discuss the thermally excited defects that terminate the plateau 
either end, among them an exotic string defect which restores three dimensionality. 



I. INTRODUCTION 

The discovery of the spin ice compounds H02T12O7 1 
and Dy2Ti2C>7 2 is one of the more remarkable events in 
the study of frustrated magnetism in the last decade. 
The name spin ice advertises their statistical mechan- 
ics at low temperature, which can - approximately - be 
mapped onto that of an Ising antiferromagnet on the 
pyrochlore lattice, which in turn is equivalent to cubic 
ice. The initial discovery stemmed from the observation 
that the large spins Jho = 8 in Ho2Ti2C>7 failed to or- 
der at any temperature despite a ferromagnetic Curie 
constant. 1 This was understood to result from the in- 
terplay of strong easy-axis single-ion anisotropy and the 
geometry of the pyrochlore lattice, which together effec- 
tively turn the ferromagnetic interaction into an antifer- 
romagnetic exchange between Ising pseudospins - which 
describe whether the moment on a given site is oriented 
inwards or outwards along the local easy axis passing 
through the site and the neighboring tetrahedra (see 
Figs. 1 and 2). 4 Later, it was pointed out that the ef- 
fective nearest-neighbor ferromagnetic exchange was in 
large part due to the effect of dipolar interactions pro- 
jected onto the manifold of Ising states. 5,6 

The antiferromagnetic interaction between the Ising 
pseudospins generates an "ice rule" - a minimum en- 
ergy configuration must involve two up and two down 
pseudospins on each tetrahedron. This ice rule does 
not, exactly as its cousin the Bernal- Fowler ice rule does 
not in the case of crystalline water, determine a unique 
ground state. Rather, there remains a residual extensive 
zero point entropy, which has been experimentally ob- 
served in the case of the Dysprosium (Jo y = 15/2 with 
ffM-B^Dy ~ 10/is) spin ice compound in good agreement 
with calculations and measurements of the entropy of 
ice. 2 Spin ice therefore offers a laboratory for studying 
the properties of water ice by proxy, but its properties 
are, of course, worth studying in their own right. For a 
review of this burgeoning field, see Ref. 7. 

It was realized early on that a magnetic field provides a 
versatile probe of spin ice, as an external field couples to 




FIG. 1: A single tetrahedron inscribed in a cube. In the py- 
rochlore lattice, the spins reside on the corners of the tetra- 
hedra. In spin ice, they are constrained to point along the 
body diagonals, d K , indicated by the short-dashed lines. The 
body diagonals define the (111) directions, the cube edges the 
(100) directions, and the bonds the (110) directions. 



the actual spin magnetic moments and thus acts on the 
Ising pseudospins in non-trivial ways. 8,9 The phenomena 
predicted here include plateaux in the magnetization 10 ' 11 
and a liquid-gas transition for fields of different strengths 
and orientations. 10 

Experimentally, the usefulness of magnetic fields was 
initially limited by the absence of single crystals, so that 
the behavior in a magnetic field had to be interpreted 
in terms of an average over all possible relative angles of 
fields and crystallites. With the advent of single crys- 
tals, this shortcoming is being removed. 12-16 In recent 
experiments on Dy2Ti2C>7 16,17 it was demonstrated that 
applying a field in a [111] direction does indeed lead to 
the predicted pair of magnetization plateaux — a low field 
plateau which retains an extensive zero temperature en- 
tropy albeit one reduced from the zero field value, and 
a second plateau at higher fields where the entropy van- 
ishes and the magnetization is saturated upon violation 
of the ice rule. 10,11 

In this paper we mostly provide a theory of the prop- 
erties of the low field [111] plateau in the T — * limit 
with some additional considerations on finite tempera- 
ture corrections and the crossovers out of the plateau 
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feet, which restores three dimensionality at low fields and 
analyze their impact on the physics at low temperatures. 

In the balance of the paper we will provide details of 
these assertions. We begin in Section II by recapitulat- 
ing the justification for using the nearest neighbor model 
and the ice rule and how they give rise to the plateaux 
of interest upon addition of a field in the [111] direction. 
We turn next to the thermodynamics and statics (Sec- 
tion III) and dynamics (Section IV) of the plateau. We 
then discuss the complications produced by the freezing 
that takes place at low temperatures, in particular with 
respect to entropy measurements (Section V), and then 
to the impact of thermally excited defects and the longer 
ranged dipolar physics in Section VI. We conclude with 
a summary. 



II. THE MODEL 



FIG. 2: The pyrochlore lattice. The cube's axes are the same 
as those in Fig. 1. 



at low and high fields. We do so within the nearest 
neighbor antifcrromagnetic model of spin ice wherein the 
low temperature limit serves to enforce the ice rule upon 
the allowed states. Further, the presence of the mag- 
netic field effects a dimensional reduction in the same 
limit — the fluctuating degrees of freedom are forced to 
live on decoupled planar subsets of the parent three di- 
mensional pyrochlore lattice which have the connectivity 
of the kagome lattice. Via a mapping derived previously 
by us in a study of frustrated Ising models in magnetic 
fields, 18,19 the remaining planar problem maps onto a 
hardcore dimcr model on the hexagonal lattice. This al- 
lows a calculation of the equal time correlations 20 - which 
are two dimensionally critical - and of the (reduced) en- 
tropy of this region which agrees well with the exper- 
iment. We next consider tilting the field weakly away 
from the [111] direction, and find that the system re- 
mains in an extended critical phase with a continuously 
drifting wavevector, 20 until it finally undergoes a con- 
tinuous phase transition, known as the Kasteleyn transi- 
tion in the dimer literature, 21 where the entropy vanishes. 
This transition has a number of interesting features, in- 
cluding the absence of any symmetry breaking, a mixed 
first/second order nature and anisotropic critical expo- 
nents. The dimer model has a height representation, and 
as discussed by Henley, 22 this leads to a natural Langevin 
dynamics for the coarse grained heights. We use this 
to write down expressions for the dynamic spin correla- 
tions in the plateau, which exhibit a dynamical exponent 
Zd = 2, although testing them is likely to be complicated 
by equilibration problems that do not affect the thermo- 
dynamics and statics. Finally, we identify the excitations 
out of the ground state manifold, which are a planar zero 
dimensional object whose condensation leads to the high 
field saturated plateau and an unusual infinite string dc- 



It is not immediately apparent that the spin ice com- 
pounds will exhibit a macroscopic low temperature en- 
tropy in zero field, let alone in a field. Indeed, a suffi- 
ciently general microscopic model for the spin ice com- 
pounds involves exchange couplings, dipolar interactions 
and a strong easy axis anisotropy. These can be encap- 
sulated in the classical Hamiltonian for unit-length spins 
S,;: 5 ' 6 



(ij) (ij) 
2 



Si ■ Sj — 3(Si ■ Tij)(Sj ■ Yij) 



+ EJ2 (d„ ( i) -Si) -gfigJ^B-Si 



(2.1) 



Here are the exchange constants and while the sum 
on (ij) runs over all pairs of sites, only a few are expected 
to be significant. The second term is the dipolar inter- 
action of strength D, where r»j is the vector separation 
of two spins measured in units of the nearest neighbor 
distance, and = Tij/\vij\. The third term is the easy 
axis anisotropy of strength E < 0, whose large magni- 
tude is crucial in these compounds and will be taken to 
infinity for the purposes of this paper, thus constraining 
the spins to point along their respective easy axes which 
we have specified by the unit vectors d K u) at site i. The 
unit cell of the pyrochlore lattice has four sites, which 
can be taken to belong to a tetrahedron of one of two 
orientations, and hence k runs from to 3 for the four 
easy axes that point from the center of the tetrahedron 
to the corner on which the site is located (see Fig. I). 
These are the also (111) directions of the underlying fee 
lattice. In the final term, we have allowed for a magnetic 
field of strength B and g/iB J is the magnetic dipole mo- 
ment of the spins. In the following, we consider fields 
along (or close to) the [111] direction. This is a threefold 
symmetry axis of the pyrochlore lattice: a field along the 
[111] direction singles out the spin (k = 0) with an easy 
[111] axis but leaves intact the symmetry between the 
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other spins (labeled k = 1,2, 3) with easy axes along the 
remaining (111) directions. 

The main difficulty in fixing the parameters in Eq. 2.1 
is lack of knowledge of the supcrexchangc, while the value 
of D can be essentially fixed via a crystal field calculation. 
It turns out that in the spin ice compounds, the effective 
nearest neighbor exchange is ferromagnetic by virtue of 
the dipolar interaction, with the weaker superexchange 
possibly being antiferromagnetic and thereby canceling 
off part of the dipolar interaction. Very little is known 
about further-neighbor superexchange, although there 
again appears to be a cancellation effect against the dipo- 
lar interactions. 5,6 

Despite these uncertainties, both experiment and the- 
ory indicate that a remarkable simplification takes place 
at moderate temperatures. If we define the pseudospins 
Ui = ±1 by whether a spin points into or out of a tetra- 
hedron on a given sublattice, i.e. we write the spins as 
Sj = crd K ^ then the accessible low energy states of 
Eq. 2.1 in zero field (B = 0) are largely governed by the 
ice rule, which requires that | J2 K a K \ = for each tetra- 
hedron. While Eq. 2.1 is believed to lead to a unique (up 
to symmetries) ground state at T = in zero field, 5,6 
this state has in fact not been observed experimentally. 
Provided Eq. 2.1 is an appropriate description, it thus 
appears that this state is dynamically inaccessible and 
irrelevant to the observed physics. 23 

The net result then is that the accessible behavior is 
captured by the greatly simplified nearest neighbor Ising 
pseudospin Hamiltonian, 

U = Jeff ^2 a i a i ~ -9Mb J ^2 B • d K(i) <Ti , (2.2) 

<ij> i 

with an antiferromagnetic J c ff. 

The ground states of this Hamiltonian for B = 
are, of course, those configurations in which ^ re a K = 
(i.e. two spins point in and two out) for each tetra- 
hedron separately. The number of these states is not 
known exactly but an estimate due to Pauling gives 
Sp/ks = (1/2) log(3/2) for the ground state entropy per 
spin which, as mentioned in the Introduction, agrees well 
with the experimental determination of the residual en- 
tropy thus providing support for the simplification. 



A. Effect of magnetic field 

The effect of switching on a field is strongly depen- 
dent on the direction of B, as first discussed in Ref. 10 
and is clear from Eq. 2.2. For instance at zero tempera- 
ture, an infinitesimal field along the [100] direction com- 
pletely lifts the degeneracy of the ensemble of spin ice 
ground states while one in the [110] direction leaves a 
non-extensive degeneracy. 

A field in the [111] direction, which is our subject in 
this paper, orders one sublattice immediately but still 
leaves a macroscopically degenerate set of ground states 



for a finite range of its values, thus producing a mag- 
netization plateau with a residual zero temperature en- 
tropy within the ice rule manifold. At a still higher field 
(gJ/iBB = 6 J e ff) the system abandons the ice rule and 
chooses the unique configuration that saturates the mag- 
netic moment in the [111] direction and thus exhibits a 
second magnetization plateau but now with no residual 
entropy. 

To see how this comes about, first note that the projec- 
tion of the total spin of a tetrahedron onto the magnetic 
field is maximized in the case of a = —1 for k = and 
a = 1 for the others. Hence at sufficiently large fields 
the system will choose the unique configuration in which 
this arrangement holds for all tetrahedra. This leads, 
however, to IX^k^kI = 2 on all tetrahedra and is thus 
in conflict with the ice constraint IX K cr K | = 0, so that 
the low field solution must be different. Instead in that 
limit one chooses a K = 1 for all the spins on sublattice 
k = as their projection onto the external field is maxi- 
mal but as the other spins have an equal projection onto 
the field, one can choose any one of these to be the sec- 
ond spin with a = — 1 needed to respect the ice rule. The 
transition between these two regimes can be located by 
computing the energies of the two arrangements. 

The pyrochlore lattice can be thought of an alternate 
stacking of kagomc and triangular planes, with the tri- 
angular planes containing all the spins of one of the four 
spin sublattices - in this case, the triangular planes of 
the k = sublattice are fully polarized and inert. Conse- 
quently, the remaining degrees of freedom live on the de- 
coupled kagome planes. Each triangle of a given kagomc 
plane has two spins with a positive projection (cr = 1) and 
one with a negative projection (cr = — 1) onto the external 
field. Such configurations are equivalent to the ground 
states of an antiferromagnetic Ising model (cr = ±1) with 
an exchange in excess of the external field (cr = { — 1,1,1} 
favored over a = { — 1,-1,1} in each triangle) 24 and as 
we show in the next section by explicit enumeration, they 
are macroscopic in number. 

While we have deduced the low field plateau (hence- 
forth simply plateau when no confusion is engendered) 
and its termination by the saturated state from the near- 
est neighbor model, its existence in experiments is further 
strong evidence for the applicability of the model and can 
be used to deduce the energy scale for the ice rule. 25 

In the next two sections we will analyze the statics, 
thermodynamics and dynamics of the plateau at low 
temperatures within the manifold of kagome configura- 
tions identified above. In Sect. VI, we will discuss semi- 
quantitatively the consequences of the inclusion of ther- 
mally excited defects that either violate the ice rule or 
are not confined to the kagome planes. We also com- 
ment briefly there on what might be missed in passing 
from Eq. 2.1 to Eq. 2.2 in our problem. 

Even with our simplifications we are left with a non- 
trivial statistical and dynamical problem that needs to 
be solved in order to compute the physical properties of 
the plateau and we now turn to this task. 
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FIG. 3: Mapping of pseudospins a — ±1 on the kagome lat- 
tice onto hardcore dimers on the (dashed) hexagonal lattice. 
Shown is the configuration favored by a field tilted slightly 
away from the [111] direction. The basis vectors used for the 
kagome lattice are also shown. 



corresponds to a tetrahedron, and hence two sites, of the 
pyrochlorc lattice, so that the entropy per spin equals 

S w 0.08077A: B , (3.1) 

which is, of course, also the value obtained in Ref. 19. 

In Ref. 16, the value obtained was 0.096 ± 0.012/cs 
per Dysprosium atom. While this work was in progress, 
another measurement has appeared, with a value of 
0.078fc_B. 17 Our value is just outside the error bars of 
the former. The fact that the former is too high suggests 
that some configurations breaking the ice rule play a role. 
Had it been too high, the implication would have been 
that a certain degree of (possibly short-range) order, pre- 
sumably due to long-range interactions, had already set 
in. If the latter, however, should turn out to be the cor- 
rect value in the end, this would be an agreement almost 
too good to have been hoped for. 

By comparison, the zero field result of Ref. 2 is 
Sq w 0.20/cs. This compares to the Pauling estimate of 
S p /k B = (1/2) log(3/2) w 0.202733 or the exact value for 
two-dimensional spin ice (for which the Pauling estimate 
is the same) oiS Lieb /k B = (3/4) log(4/3) « 0.215762, so 
that the decrease due to the applied field is by a factor 
of 2.5 - 2.7. 



III. PLATEAU: THERMODYNAMICS AND 
STATICS 

In the last section we noted that the allowed spin con- 
figurations in a single kagome layer in the plateau are 
equivalent to the ground states of the Ising antiferromag- 
net on the kagome lattice. We have previously consid- 
ered this problem and shown that the ground states are 
in correspondence with the configurations of the exactly 
soluble problem of the dimer model on the honeycomb 
lattice, 18 a mapping rediscovered by Udagawa et al. 19 
The triangles of the kagome lattice form a dual hexag- 
onal (honeycomb) lattice, whose bonds are the sites of 
the kagome lattice. For each spin with positive projec- 
tion onto the field, color in the corresponding link of the 
hexagonal lattice. As each triangle has exactly one such 
spin, each site of the hexagonal lattice has exactly one 
colored link emanating from it. By calling the colored 
link a dimer, one thus establishes an exact one-to-one 
correspondence between the configurations of a hardcore 
dimer model on the hexagonal lattice and the spin ice 
states in a weak [111] field. 



A. Entropy 

The entropy of the dimer model on the hexagonal lat- 
tice is well known, having been first computed as the 
entropy of the equivalent triangular lattice Ising anti- 
ferromagnet at T = 0. The latter has an entropy of 
Sa = 0.32306fcs per site. This corresponds to an entropy 
of So = S&/2 per site of the dimer model. Each triangle 



B. Correlations 

The dimer model describing the plateau has a range of 
further interesting features in addition to its nonvanish- 
ing zero point entropy. Most strikingly, its correlations 
are critical, decaying as 1/r 2 at large distances, r. 

In detail, consider the connected pseudospin correla- 
tion function 

c«x(r) = (a A (rK(0)) - (a x ) (a K ) , (3.2) 

where r labels the location of the tetrahedron and the 
Greek letters the location of a pseudospin in the tetrahe- 
dron. This is simply related to the correlation functions 
of the real spins, C kX {t) = (S A (r)S K (0)) - <S A ) (S K ). For 
instance, for the components of S along the [111] direc- 
tion, 

= (S/3) 2 (-3) s «-°+^°c kX . (3.3) 

where the factors of 3 are due to the different projections 
of the inequivalent easy axes onto the [111] direction. 
Similarly, the full spin-spin correlation function is given 
by 

C kA = -(S/V3) 2 (-3) 5 - a c«a. (3.4) 

In the plateau region, ao = (ctq) = — 1 everywhere, 
so that Co, a = 0. The nontrivial correlations involve 
only k > 0, that is to say spins in the same kagome 
planes. These correlations can be calculated following 
Ref. 20. We have tabulated the short distance correla- 
tions in Fig. 4. 
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FIG. 4: Short distance correlations, 10 x ci«/4 of the 
(pseudo)spin in the bottom left hand corner, marked by a 
solid dimer. Positive correlations are indicated by dashed 
dimers. This plot uses the same normalization conventions 
as that of Tab. I in Ref. 20, hence the factor of 1/4; in the 
convention of the dimer model, the correlation at the origin is 
— 1/9 — »• —11111. Recall that the dimers occupy the links of 
the hexagonal lattice, the midpoints of which are the kagomc 
lattice sites. 

The correlations decay algebraically at long distances. 
The two independent correlators are, 

c n( r ) ~ TTT^ [cos(47nc/3) - cos(20)] (3.5) 
ci 2 (r) ~ -4-7 [cos(4ttx/3 + 4tt/3) - cos(26> + 4rr/3)] . 

Here, r is a Euclidean coordinate vector for the kagomc 
lattice, with r = |r| being the distance between two tri- 
angles of the kagome lattice, and tan# = y/x, see Fig. 3. 
This asymptotic behavior, involving a sum of oscillations 
at wavevector q x = 4n/3 and a dipolar piece, can be 
readily obtained by means of the height representation 
formulae listed in the next section as well. 

As a consequence of the first term in brackets in 
Eq. 3.6, one would therefore expect a peak in the 
Fourier transform of the structure factor at wavevector 
±(47r/3, 0). Here we have used the lattice constant, twice 
the pyrochlore nearest neighbor distance, as the unit of 
length (see Fig. 3). 

The corresponding peaks at the four symmetry related 
locations are obtained by the appropriate addition of re- 
ciprocal lattice vectors, 2tt(1, — 1/\/3) and 2n(0,2/>/3). 
Note in particular that 27r(2,0) is the reciprocal lattice 
vector relating the peaks at — (47r/3)x and (87r/3)x. How- 
ever, in Fig. 5, the peak at the latter location is absent. 
This happens because the 'form factor' of the unit cell 
has a zero at (8n/3)x, as can be verified directly from 
Eq. 4.4. In Fig. 6, this effect is reversed in that the peak 
at (87r/3)x is the stronger one; the peak at (4*7r/3)x, al- 
though present, is not visible on the contour plot for the 
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FIG. 5: The Fourier transform of the pseudospin correlations, 
c, in the kagome planes, obtained from a finite system con- 
taining 9604 sites. q x ,q y range from — 4ir to 4tt. In addition, 
there is a peak at q = and the reciprocal lattice vectors due 
to the finite average moment induced by the field. Note the 
logarithmic peak at (47r/3, 0) and the symmetry related po- 
sitions. Together, they should describe the differential cross 
section found in polarized neutron scattering with the neutron 
spin pointing along the [111] direction. Light regions denote 
strong scattering. 



system size considered as it is almost an order of magni- 
tude weaker. 

These are not true Bragg peaks, as there is no long 
range order. Indeed, as the power law decay of the pseu- 
dospin correlations is rather rapid, r~ 2 , their intensity 
grows only logarithmically with the planar system size. 
Similarly, the intensity decreases logarithmically as one 
moves away from the center of the peak. The second 
term, although of equal amplitude, does not lead to a 
feature with macroscopic intensity, as no finite fraction 
of its weight is concentrated on any one wavevector. In 
Fig. 5 we plot the absolute value of the Fourier trans- 
form of the full pseudospin correlation function, which 
exhibits these features and is detectable by polarized neu- 
tron scattering. In Fig. 6 we plot the cross-section for 
unpolarized neutrons; the difference in the two figures re- 
flects the non-trivial relation between the spins and the 
pseudospins. Both figures omit the magnetic Bragg peaks 
that will arise from the static magnetization produced by 
the applied field, and are obtained for zero out-of-planc 
wavevector transfer. 
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E? = -(^ b BJ/3)cti 



FIG. 6: The Fourier transform of the correlation of the spins 
components perpendicular to the in-plane wavevector. Details 
as in the previous figure. The quantity plotted here is also the 
differential neutron scattering cross section for unpolarized 
neutrons. 



C. Kasteleyn transition in a tilted field 

A broader view of the critical correlations in the honey- 
comb dimer model is obtained by generalizing it to allow 
for unequal fugacities for dimers of different orientations. 
As shown by Kasteleyn, 21 the equal fugacity point sits in 
a critical phase which borders a "frozen" phase with van- 
ishing entropy that is reached by an unusual transition 
that bears his name. If zi, z 2 and Z3 are the fugacities of 
the three sets of dimers, the transition takes place when 
the fugacity of one set equals the sum of the other two, 
say z\ = z 2 + z$. For z\ > z 2 + z 3 a unique configura- 
tion survives (shown in Figure 1). It is interesting to ask 
whether this phase transition can be realized in the spin 
ice problem. It turns out that this can be done rather 
simply by tilting the field. 

To see this consider tilting the applied field away from 
the [111] direction so that it acquires an enhanced com- 
ponent in the [-1-11] direction, which is the easy axis 
of sublattice k = 1: B = B (cos<^)[lll]/v / 3 + sincb [-1- 
12] /VfJ), so that the angle the field makes with the [111] 
direction is given by <j>. This keeps the other two of the 
three kagome spin sublattices (/c = 2, 3) equivalent and 
singles out the k = 1 sublattice. To leading order in the 
tilt angle, spins on sublattice k = do not experience a 
change in energy, whereas spins on the other sublattices 
do: 



COS ( 



E 2,3 



-(gfi B BJ/3)a 2 ,3 



cos 



- 2\/2sin(? 
4> + V^sini 



(3.6) 



As the dimer fugacities are z K = exp[2E^ /(fcgT)], it fol- 
lows that the effect of the tilted field is to make them 
unequal - specifically, to privilege the occupation of ver- 
tical dimers over the other two orientations in Fig. 3. At 
zero temperature z± is infinitely bigger than z 2 or Z3 at 
any tilt angle and the system is deep in the frozen phase, 
which is to say the energy gain is all there is and we ob- 
tain just the so-called staggered configuration shown in 
Fig. 3. 

At nonzero temperatures, or finite fugacities, however, 
the gain in energy must compete with the loss of en- 
tropy, both extensive, to effect a gain in free energy and 
we obtain a finite range of stability for the critical phase 
terminated by the Kasteleyn transition. From the cri- 
terion z\ = z 2 + z 3 we can deduce a critical tilt angle 
4>c set by ksT = (2>/2/ln2)g fisB J sin</>, at which the 
transition occurs. Note that the transition temperature 
is proportional to the in-plane field strength, Bsm<fi, so 
that the experiment can, in principle, be done at T <C B 
and when the tilt angle is sufficiently small to justify our 
neglect of 0{<p 2 ) terms. In the following, we express the 
dependence on the various parameters via z = z 2 jz\. 

Various predictions follow from this analysis: 
(a) The Kasteleyn transition involves a critical vanishing 
of the entropy 



S ~ (0 C - 0) 1 / 2 



(3.7) 



that can be detected via standard thermodynamic mea- 
surements. In equilibrium this implies a significant sig- 
nature in the tilt specific heat, C, in the form of a diver- 
gence, 



(0c - <£) 



-1/2 



(3.8) 



but freezing is likely to complicate such a direct measure- 
ment as we discuss in Section V. 

(b) The expectation values of the Ising spins for z > 1/2 
is given by 



(01) = — H arcsin[yi 

<<72> = <*3> = (l-<*l»/2. 



4z 2] 



(3.9) 



The magnetization in the [-1-12] direction, m±, being 
proportional to {<T\), it follows that it deviates in the 
critical region from its saturation value, m^, as 



- 



{4>c - <P) 



(3.10) 



EE 



g[iBBJ<Jo cos ( 



This expression holds to the left of the critical point (z > 
1/2, see Fig. 7). To the right, there are no fluctuations, 
and (<7 2 ) = (cr 3 ) = - (01) = 1 

The correlations remain critical but change continu- 
ously as B is tilted. For example, the equation for the 
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same sublattice connected correlations, Eq. 3.6, is gener- 
alized to 20 ' 26 

cn(r') = -^ti [«»(2a:/e e ) - cos(20')] ■ (3.11) 
2-7TT 

Here, r' 2 = x 2 + V. with 

l/£ x = 2arcsin A /l - l/4z 2 (3.12) 
1/^ = (4z/V3)v/l - l/4z 2 arcsi nv /l - l/4z 2 , 

with z = z 2 /zi = exp[2(£?f - E^)/(k B T)} and tand' = 
i x / £x) / (y / £ y ) ■ From this we observe that: 

(c) The location of the peak in the structure factor, which 
remains logarithmic, is given by ±(2/£ x ,0), so that it 
drifts continuously from (Air/3)x to the center of the Bril- 
louin zone, which it reaches at the phase transition. Ob- 
servation of this drift with field tilt should be a good flag 
of the unusual critical phase. 

(d) The scattering pattern is reduced in symmetry - 
the applied field reduces the six-fold rotational symme- 
try of the lattice to a two-fold one. In particular, this 
leads to anisotropic scaling at the Kasteleyn transition in 
which there are two diverging correlations lengths along 
(£x ~ (0 ~ c 0) -1 ^ 2 ) and transverse (£ y ~ ((f) c — (j))^ 1 ) to 
the in-plane field, whose ratio £, y /£, x also diverges as one 
approaches the transition, z — > l/2 + . 

(e) Finally we note that the transition is asymmetric. 
On the side z — > 1/2 - , no fluctuations are present, so 
that the transition has an asymmetric first/second or- 
der appearance. However, the latter property is strictly 
dependent on the hardcore condition on the dimers and 
tetrahedra violating the ice rule will allow some fluctua- 
tions even beyond the transition, see Sect. VI. 

IV. PLATEAU: DYNAMICS 

We now turn to the dynamical correlations in the 
plateau continuing to assume that the system explores 
only its ground state manifold; we will return to the va- 
lidity of this approximation in Section V. Prima facie, 
finding the time dependent correlations seems a difficult 
task since the configurations are characterized by a lo- 
cal constraint, which we have compactly represented by 
the hard core dimer mapping. Nevertheless, this can be 
done at long wavelengths and low frequencies, following 
the ideas of Henley on the dynamical correlations of crit- 
ical dimer models, 22 which we apply to the honeycomb 
lattice in the following. Henley's basic insight is that the 
dimer configurations on bipartite lattices have a height 
representation whose fluctuations are unconstrained at 
long wavelengths. For the statics this has been known 
since the work of Ref. 27 (see also Refs. 22,28 for a con- 
cise introduction) and the extension to dynamics leads 
naturally to a Langevin dynamics for the heights. The 
resulting theory is Gaussian and exhibits dynamic scaling 
with the dynamic exponent Zd = 2. We now give brief 
details of this analysis. 




-0 .25 



FIG. 7: Magnetization of the spins in the kagome planes in 
the [-1-12] direction (thin line) and inverse correlation lengths 
(thick lines) in the x direction and y direction (in black). The 
former is normalized with respect to the saturation magne- 
tization for B > B c , m sat = (4V2/3)<7^s J ■ Saturation for 
B — > — oo is half this value (and negative). The inverse corre- 
lation lengths are normalized to their zero field value of 2n/3. 
Note that they vanish with different powers at the transition. 
The x coordinate is given by (2^/2/ \n2)gfi B B J sin (j)/(k B T), 
so that the critical point is located at 1. 

First, we provide a description of the relevant height 
model. Microscopically, this involves a map between 
dimer configurations and the configurations of a surface 
specified by giving its local height above the dimer plane. 
The microscopic heights are a set of integers, defined on 
the sites of the triangular lattice dual to the hexagonal 
lattice the dimers reside on. The height changes by +2 
(—2) if one crosses a dimer when going from one site to 
its nearest clockwise neighbor on an up- (down-) trian- 
gle. If no dimer is crossed, the change is —1 (+1). This 
provides a mapping of dimers onto heights. The dimer 
density, rid, is thus given by nd = (V lat /i + l)/3, where 
yiat denotes the lattice derivative corresponding to the 
rules defined in this paragraph. 

In the coarse grained, continuum theory, this micro- 
scopic expression indicates the identification, 

n d =~@>V)h+~, (4.1) 

where e is a unit vector perpendicular to the orientation 
of the dimer. This is however, not the full expression, 
even at leading order. Upon coarse-graining, a second 
non-trivial term appears in the expression for nd, which 
reflects the important fluctuations near the characteristic 
wavevector of the flat states - this is the analog of the 
staggered "2fc^-" piece that appears in the bosonization 
of one dimensional quantum fermion systems. This piece 
can be identified by noting that the mapping of dimers 
onto heights is one-to-many: a shift of the height by 3 
units returns the same dimer configuration, and thus the 
operator must be invariant under this operation. 2 One 
thus obtains for the dimer densities n K : 
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1 1 

m-- = -d x h + (exp(2nih/3)exp(4:irix/3) (4.2) 

n 2 -^ = ^(-^a x + ^9j,)ft + Cexp(27ri/i/3)exp(47ria;/3 + 47ri/3) 

n 3 - i = ^(-i^ - ^a s/ )/i + Ccxp(27r«/i/3)cxp(47ria;/3-47r«/3). 

where the normalization £ = l/(2na) involves a short distance cut off a. There are, of course, corrections from less 
relevant operators which we have not considered here. 

To calculate the static dimer correlators, one uses the fact that the heights fluctuate in a Gaussian manner in 
equilibrium, 



H 



J d 2 r^\Vh\ 2 (4.3) 



(K = tt/9 for the honeycomb lattice), whence the height correlator is given as (h(r)h(O)) = — ln(r/a)/(2nK). From 
these we find the asymptotic correlations, 



c v ( r ) = ^2 [cos(4^/3 + M? - 0/3) - cos(2# + + j - 2)/3)] , (4.4) 



in agreement with Eqs. 3.6. One sees that the two pieces 
in the dimer correlators arise from the "uniform" and 
"staggered" pieces of the representations given above. It 
is also straightforward to check that the structure factor, 
at this level of approximation, gets no contribution from 
the uniform pieces and consists entirely of the logarithmic 
peaks at ±4tt/3x and related points. In addition, the 
extinction of the peaks at ±87r/3a; in Fig. 5 also follows 
from Eq. 4.4. 

To obtain the dynamical correlations, we note that 
the long wavelength, low frequency dynamics for a 
generic local dimer dynamics will be governed by Hen- 
ley's Langevin equation 22 



dh(r) 
dt 



= -r- 



SH 
Sh(r) 



CM) 



(4.5) 



where T is a kinetic coefficient set by microscopies and 
the noise £(r, t) obeys 



(C(r,t)C(iV)>=2r5(r-r') ( J(t-f) . 



(4.6) 



As this is again a Gaussian theory, it follows that the 
only non-trivial correlator of the heights is the two-point 
function, 



h q {t)h- q (0) 



1 



Kq' 



exp[-X(q)t] , 



(4.7) 



where h q (t) is the height configuration at wavevector 
q = (q x ,q y ) and time t. The relaxation rate for the 
modes with wavevectors of magnitude q is given by 
X(q) = TKq 2 , which implies a critical dynamics with 
z d = 2. 

The dynamic correlations can now be obtained from 
this expression in the same manner as the static one. 



For example, the uniform piece of the same sublattice 
correlator equals 

KW«i-.(0)) u = jLexp[-A(g)i], (4.8) 
which yields the further Fourier transform, 



Kq 2 q^ + w 2 ' 



(4.9) 



As in the case of static correlations, the structure factor 
gets no contribution from such uniform pieces. 

The non-zero contribution then comes from the stag- 
gered piece which is first calculated in real space as the 
vertex operator correlator, 

(ni(r,i)ni(0,0)) a = ^{ e ^ ix / 3 ( e 2vih ^/ 3 e- 2vih ^/ 3 ) + c.c.} 



„ a9 / AttxX ( 47T 2 
2C 2 cos^— joxp^- — 



C(r,t)\ , 



where 

C(r,t) - ([h(r, t) - h(0,0)} 2 ) /2 = (4.10) 
JJ(0^ [1 - CXp( - r ^ )COS(q ' r)] - 



In the scaling limit, (r, t) — ► oo with r 2 /t fixed, this can 
be written in the scaling form, 



1 



/ 4:7TX\ 



( ni (r,t) ni (0,0)). = ^cos y—jg y—j 

(4.11) 

where the scaling function is given in terms of the incom- 
plete Gamma function as 



l(x) = e 



-r(o,x/4) 



(4.12) 
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and exhibits the asymptotics 

< \ / e7x / 4 X <C 1 (A1*\ 

9[X) ~ \4exp(-x/4)/a; a; » 1 ' 1 j 

The former encodes the autocorrelation 

^1(0,4^1(0,0)). = ^^, (4.14) 

where 7 = 0.5772... is the Euler-Mascheroni constant. 

The remaining task is to obtain the Fourier transform 
of Eq. 4.11 which does not appear possible in closed form 
and will therefore probably to be accomplished numer- 
ically if desired. However, the essential features can be 
deduced as follows. 

First, the Fourier transform will still be peaked about 
±4-7r/3x and symmetry related points. Second, if we mea- 
sure momenta from each of these values, the result ex- 
hibits the scaling form 




(4.15) 



Third, one can show that 

g(x) ~ ( V 2 *<<J (4.16) 
y v ' \ c/x x » 1 y J 

with some constant c and that the corrections about 
cither limit arc analytic. Together, the lasts two fea- 
tures imply that fixed frequency cuts will exhibit peaks of 
height (2|u;|) _1 , finite with divergent system size, whose 
widths will exhibit the characteristic Zd = 2 scaling, 
Aq <~ ^/lu/TK. The complimentary fixed q cuts will 
exhibit a diffusive peak at uj — of height c/ (TKq 2 ) and 
width Auj ~ TKq 2 . 

It is worth noting that in taking the scaling limit we 
have kept all information relevant to long wavelengths 
and low frequencies but if we attempt to reconstruct the 
equal time correlator we will find a spurious ultraviolet 
singularity. Likewise the large frequency behavior at a 
fixed q will be softer than the l/\co\ dependence implied 
by the scaling form. 

V. FREEZING 

This is a good place to note an important subtlety in 
making contact between our analysis, and indeed all the- 
oretical work on ice and spin ice, and the experimental 
systems. This is the feature that both ice and spin ice 
exhibit diverging relaxation times (set by the tempera- 
ture dependent T in our formalism) at low temperatures 
which overtake the timescale of experiments so ergodicity 
is lost. For spin ice the evidence for this comes from the 
experiments of Refs. 29,30, which report a strong slow- 
down of the dynamics setting in around 1-2 K, a signa- 
ture of which is the appearance of hysteresis in magneti- 
zation measurements. Consequently we need to examine 



whether the equilibrium computations of the this paper 
represent measurable quantities. 

The good news is that the thermodynamic and static 
quantities are indeed still measurable. For the magne- 
tization and the static structure factor this is a con- 
sequence of self-averaging in the sample - with proba- 
bility one these quantities are the same for a configu- 
ration picked at random as they are for the entire en- 
semble of ground states. This in turn comes from two 
sources. First, even in a frozen three dimensional configu- 
ration, the different Kagome planes effectively give differ- 
ent members of the equilibrium two dimensional ensem- 
ble. Second, even in a given plane we get self-averaging. 
For example, the spin-spin correlation function at a fixed 
separation, averaged over the location of the spins in a 
configuration picked at random, converges to its ensem- 
ble averaged value in the limit of infinite system size; the 
algebraic correlations in our problem lead to at best a 
(log N) 1 / 2 correction to the 1/VN dependence expected 
for the fluctuations in a system with iV sites. As the 
structure factor involves exactly this average, all is well 
on that front. The same holds for the magnetization, 
measured as the moment frozen into a field cooled sam- 
ple. 

The story with the entropy is different. Indeed it is 
worth emphasizing the remarkable fact that experiments 
measure an entropy associated with a macroscopic de- 
generacy of ground states even as the system settles into 
just one of them (or a sub-macroscopic number since lo- 
cal fluctuations presumably do survive even as large scale 
rearrangements are frozen out). The contradiction with 
the statistical mechanical view of entropy as the loga- 
rithmic volume of phase space explored is resolved when 
one notes that the experimental determination consists 
of starting with the known entropy of the paramagnetic 
high temperature state and integrating down with the 
measured heat capacity. At issue then is whether the 
freezing substantially affects the ratio of heat capacity 
to temperature over the temperature range where it is 
significant. For the ice problems, the spectrum involves 
a finite gap to making a defect above the ground state 
manifold. Consequently, at temperatures below this gap, 
which is also where freezing takes place, the heat capac- 
ity is exponentially small in the temperature, whence the 
freezing hardly affects the entropy determination. 31 In 
our problem this implies that field cooled measurements 
of the heat capacity will allow determination of the ther- 
modynamic entropy inclusive of tilted field values. 



VI. THERMAL AND ANALYTIC DEFECTS 

Thus far our analysis has assumed that the only acces- 
sible configurations belong to the ground state manifold 
of the pscudospin Hamiltonian. To make contact with ex- 
periments we need to examine the effects of relaxing this 
restriction. In this section we do this, thereby obtain- 
ing some insight into the low and high field boundaries 
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of the plateau and also comment on a couple of other 
salient limitations of our analysis. 

As noted earlier, at T — simple energetics shows that 
the plateau extends over < g^sJB < 6J c ff, giving way 
at zero field to the full spin ice ground state manifold and 
to the right to the fully saturated state. At finite tem- 
peratures the plateau state is no longer field independent 
but will instead evolve, especially near the transitions. 
At low temperatures we can gain insight into this evolu- 
tion by examining the thermally excited defects that will 
dress the critical dimer state that we have discussed in 
this paper. 



A. Monomer defects 

The first defect to consider increases the local magne- 
tization and it is the condensation of such defects which 
terminates the plateau at its high field end. The local 
minimum energy process to consider is one in which a 
down pseudospin in a kagome plane is converted to an 
up pseudospin so that all spins of the two triangles that 
share it are now aligned with the field. Such a process 
violates the ice rule as there are now two tetrahedra with 
^ fi iT K ^ 0, and takes us out of the ground state mani- 
fold. A single flipped spin in fact corresponds to a pair of 
defects, which is most easily seen in the dimer representa- 
tion where it corresponds to two monomers on adjacent 
sites of the hexagonal lattice. The two partners of the 
pair can be separated by moving one of the defects, on 
an 'up' triangle, say, to a neighboring up triangle. This is 
done by flipping two spins on an adjacent 'down' triangle, 
namely the a = — 1 spin and the spin it shares with the up 
triangle. This puts the original up triangle back into the 
spin ice ground state at the expense of violating the con- 
straint on the up triangle sharing the spin with the down 
triangle, ft follows then that the energy cost of flipping 
the spin is the creation energy 2£ m — 4J ff — 2gfisJB/3 
of two defects. This energy vanishes exactly at critical 
field g^bJBc = 6J e s which separates the two plateaux at 

rp _ Q 10,11 

At finite but low temperatures, the system contains 
a finite but small density of these defects whose sepa- 
ration will set a correlation length and cutoff the criti- 
cal singularities of the parent dimer state. Naively, we 
might anticipate £ 2 <~ l/n m <~ exp(£ m //cBT) but there 
is a pseudo- Coulomb (logarithmic) entropic interaction 
between them that modifies this dependence. The ex- 
act dependence can be computed by an energy-entropy 
balance argument that is equivalent to a tree level renor- 
malization group computation. 32 Consider a system of 
area A and let Z(ri, r 2 ) be the number of configurations 
of the dimers (spin background) in the presence of the 
two monomers (defects) held fixed at positions ri and 
r 2 while Z is the number of configurations of the dimers 
with no monomers present. Then the free energy cost of 



introducing two defects is 

AF = 2£ m -Tlog J d 2 n J d 2 r 2 Z(r 1 ,r 2 )/Z . (6.1) 

The ratio Z(r 1 ,r 2 )/Z can be computed by height rep- 
resentation theory by noting that monomers on the 
two sublattices correspond to a height mismatch of ±3 
when encircled. The operator identification described in 
Rcf. 27 then implies 

Z(n,r 2 )/Z~ 1/V|ri -r 2 | , (6.2) 

which is the same decay first described in Ref. 33 for the 
closely related square lattice dimer problem. With this 
in hand, it is easy to see that AF < when the system 
size £, which we now identify with the correlation length 
is given by 

f 2 ~V«™~<* p(^). (6.3) 

B. Termination of the plateau by monomers 

At a fixed location in the plateau the above formula 
will describe the asymptotic low temperature approach 
to the purely dimer manifold. At a fixed temperature 
though this analysis will break down near B c where a 
treatment of the statistical mechanics of large numbers 
of defects needs to be devised. We expect to address 
this problem in more detail elsewhere and here we will 
content ourselves with three remarks. 

First, matters simplify in a scaling limit T — > and 
B — > B c with (B — B c )/T fixed. In this limit we can ig- 
nore all spin configurations save those consisting of dimer 
configurations "doped" with some number of monomers. 
The remaining problem is the non-interacting monomer- 
dimer problem and hence the interpolation between the 
two plateaux as a function of B is a crossover and not 
a phase transition. 34 Second, at the transition field, this 
leads to an equal weight sum over all monomer-dimcr 
configurations. The entropy at this point is then higher 
than it is in the low field plateau before it turns around 
and then heads for zero deep into the high field plateau. 
Third, the transition point exhibits a temperature inde- 
pendent ensemble in this treatment which should lead to 
a crossing point for the magnetization isotherms. Above 
a critical temperature, the data 1 indeed exhibit a maxi- 
mum in the entropy and a crossing point for the magne- 
tization isotherms. Below this temperature the crossover 
appears to turn into a first order transition at which point 
the entropy plummets with temperature and the mag- 
netization develops a discontinuity. 17 Prima facie this 
appears to be a puzzle for the nearest neighbor model 
considered in this paper, although it is possible that a 
purely mean field treatment of the longer ranged pieces 
of the dipole interaction omitted here renormalize B suf- 
ficiently to turn the sharp low temperature crossover into 
a transition. 
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C. String defects 

The second type of defect to consider is responsible for 
decreasing the magnetization towards the low field end 
of the plateau. As in this limit we must preserve the ice 
rule, decreasing the magnetization requires that we flip 
a spin on the triangular sublattice k — while satisfy- 
ing the ice rule by choosing a second spin in the kagomc 
plane to have a = — 1. Interestingly, this is not enough 
since the k = spin is shared by another tetrahedron and 
so on. Indeed, one can sec quite generally that it must 
be infinite in length. This follows from the observation 
that the local ice rule leads to the global property that 
all [111] triangular planes have the same magnetization, 
which is equal and opposite to that of all kagome [111] 
planes. 35 As the magnetization of the triangular [111] lay- 
ers is saturated, reducing it by flipping one of its spins 
in one layer requires flipping one spin in all of the other 
layers at the same time. The energy of such a defect, 
£ Sl is thence most conveniently quoted per (kagome and 
triangular bi-) layer. As it involves antialigning a spin 
in the triangular and one in the kagome layer with the 
field, we have £ s — 8g^ B JB/3. Remarkably, despite the 
energy cost proportional to the linear system size, L, it 
is still entropically favored in a large system. To see this, 
note that such a defect corresponds to inserting a sur- 
plus dimer, violating the hard core condition, into each 
kagome plane, which connects a (say) up triangle above 
which a spin on sublattice n — is flipped with a down 
triangle below which the next flipped k = spin is lo- 
cated. As in the case of the pair of monomers defects, 
the pair of triangles can again be separated into two dis- 
tinct defects — in dimer language into two sites with two 
dimers each. If the separation of these sites were to cost 
no (in plane) entropy, one would be free to choose which 
of the A spins in the triangular layer to flip, thereby en- 
dowing the defect with an entropy of £2 = hi A per layer. 
For a sufficiently large system, it would therefore always 
be free energetically favorable to generate such a defect. 

The actual density of such defects is lowered by the 
same in-(kagome)plane entropic mechanism discussed for 
monomer defects. Again we appeal to height represen- 
tation theory to find that sites with two dimers carry 
charge ±3 so that the entropic interaction between them 
is the same as for two monomers. This implies that 
per layer the entropic gain from being able to pick the 

separation of the defects grows as log J A r~ x l 2 rdr <~ 
(3/4) In A. From this we deduce that a cylinder of cross 
sectional area A first nucleates this string defect when 
3/4 In A = (8/3)gfi B JB whence we expect the area den- 
sity and hence transverse correlation length set by 

£ 2 ~ l/n s ~ exp[32g(j, B JB/9k B T} (6.4) 

at low temperatures. This exponential dependence will 
then determine the approach of the magnetization to its 
plateau value at a fixed low field as temperature is low- 
ered. 



Again, the proliferation of such defects at low fields but 
fixed temperature requires a different treatment, involv- 
ing a linear response calculation about the full spin ice 
manifold, which we will discuss elsewhere. In this regime 
all relevant energies are set by the field so that physical 
quantities will be functions of B/T alone. We expect 
then that the magnetization curves will collapse with a 
finite slope at the origin when plotted as a function of 
B/T. 

We can draw one further inference from our computa- 
tion of the defect densities. By equating the activation 
energies of the two defects we can identify the field at 
which their densities cross at the lowest temperatures — 
this will also be the field at which the magnetization 
isotherm crosses the zero field value of the magnetization 
at low temperatures and hence a second crossing point. 
This yields a field g/i B JB = (18/31) J e g which is about 
a tenth of the critical field between the plateaux. 

How does the presence of such defects alter the results 
we have described above? Fundamentally, their presence 
will of course make itself known as a deviation from the 
'exact' result; in particular, the smallest of the defect 
induced finite correlation lengths will determine the cut- 
off at which, for example, the logarithmic peaks in the 
neutron scattering stop growing. 

As for the Kasteleyn transition, both types of defects 
will inevitably smear out the fluctuation-free regime and 
therefore the mixed first/second order nature of the tran- 
sition. Monomer defects can be exponentially suppressed 
by lowering the temperature (compared to J). As one 
lowers the temperature at small fields B/T <C 1, the 
angle sin (f> at which the transition takes place decreases 
inversely with B/T, whereas the density of string defects 
is exponentially suppressed. By achieving an improved 
angular resolution, the crossover from Kasteleyn behav- 
ior to a more conventional second order phase transition 
could thus be reduced. 



D. Disorder and dipoles 

Finally we turn to two significant limitations of our 
analysis in this paper. First, actual samples are likely to 
contain structural defects due simply to chemical disorder 
such as vacancies or interstitials affecting site occupancy 
or exchange paths. We are not aware of a determination 
of the density of such defects, although for Hciscnberg 
spins on the related SCGO lattice, there have been both 
experimental 36,37 and theoretical 38 attempts to deter- 
mine the density of vacancies from thermodynamic 36 ' 38 
or NMR experiments. 37 As the chemical defect density in 
single crystals tends to be higher than in powder samples, 
this might be a not insubstantial effect in this context. 

The second important feature omitted from the 
nearest-neighbor spin ice model are the effects of the 
long-range dipolar interactions beyond the nearest neigh- 
bor piece, which are sizeable due to the large spin of the 
Dysprosium ion. We have already alluded to one possi- 
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ble effect in our discussion of the transition between the 
two plateaux, namely that the polarization of the spins 
may require a self-consistent treatment of the field B that 
acts upon them. While this is always necessary when a 
macroscopic magnetization is present, in our case the is- 
sue is somewhat more delicate since the largest piece of 
the dipolar interactions has already been accounted for 
in the nearest neighbor model. 

On a fundamental level, however, the long-range dipo- 
lar interactions do not seem to lead to a significant inter- 
plane ordering effect, as this would have reduced the en- 
tropy determined in the experiment. This may, however, 
be a consequence not of the precise thermodynamic be- 
havior of the spin ice Hamiltonian in a field, but rather 
an indication of the magnet's inability to access its true 
ground state in the presence of energy barriers as dis- 
cussed in Section II. 



VII. SUMMARY 

The application of a field in the [111] direction to the 
spin ice compounds leads, by a reasonable set of approxi- 
mations, to an elegant dimensional reduction of the three 
dimensional problem onto a set of decoupled two dimen- 
sional problems. Fortunately, the resulting two dimen- 
sional problem is one of planer dimers and hence is ex- 
actly soluble, so that the statics and thermodynamics 
can be determined exactly. While the computed entropy 
has already been measured, the predictions for the cor- 
relations can be tested by scattering. Also testable are 
thermodynamic and static predictions for a Kasteleyn 
transition upon tilting the field in the [-1-11] direction 
and for the dynamic correlation in the plateau. Finally 
we have sketched a theory of the finite temperature mod- 
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